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Abstract 

A new network model is proposed to describe the 1/ f a resistance noise in disordered materials for 
a wide range of a values (0 < a < 2). More precisely, we have considered the resistance fluctuations 
of a thin resistor with granular structure in different stationary states: from nearly equilibrium 
up to far from equilibrium conditions. This system has been modelled as a network made by 
different species of resistors, distinguished by their resistances, temperature coefficients and by the 
energies associated with thermally activated processes of breaking and recovery. The correlation 
behavior of the resistance fluctuations is analyzed as a function of the temperature and applied 
current, in both the frequency and time domains. For the noise frequency exponent, the model 
provides < a < 1 at low currents, in the Ohmic regime, with a decreasing inversely with the 
temperature, and 1 < a < 2 at high currents, in the non-Ohmic regime. Since the threshold current 
associated with the onset of nonlinearity also depends on the temperature, the proposed model 
qualitatively accounts for the complicate behavior of a versus temperature and current observed 
in many experiments. Correspondingly, in the time domain, the auto-correlation function of the 
resistance fluctuations displays a variety of behaviors which are tuned by the external conditions. 

Keywords: Resistor networks, fluctuation phenomena, 1/f noise, disordered materials, non-equilibrium pro- 
cesses 
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I. INTRODUCTION 



The analysis of resistance fluctuations has proved to be a very powerful^tooljor^probing 
various condensed matter systems Q, [2I, [J, including nanostructures 



and disordered materials, like conductor-insulator composites 
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One of the most relevant features of the resistance noise 
lies in its dependence on frequency. Many condensed matter systems display the so called 
Lorentzian noise , which is characterized by a power spectral density of the resistance 
fluctuations scaling as 1/ f 2 at high frequencies and becoming flat below a corner frequency 
f c . A behavior associated in the time domain with an exponential decay of the correlations 
and thus with a well defined characteristic time, r (correlation time) [J, Q]. 

On the other hand, it is well known that many other condensed matter systems 

exhibit 1/f resistance noise, i.e. a spectral density scaling at low frequencies as 1/ f a with 
a ~ 1, thus a noise associated with a non-exponential decay of the correlations in the 
time domain. The ubiquitous presence of 1/f noise in a large variety of phenomena of 
very different nature has given rise to many attempts to explain it in terms of a universal 



law 



Moreover the link between 1/f noise and extreme value statistics has been 



also investigated [3l|, |32j. A simple way to obtain an 1/f spectrum is by superimposing a 



arge number of Lorentzian spectra with an appropriate distribution of the correlation times 



In some cases, this distribution can be derived from the distribution of some variable 



on which the correlation times themse 
Mc Whorter and Dutta et al 



yes depend, as in the case of the pioneering works of 

3J 



, [33| . In particular, these authors proposed 
that the origin of the 1/f noise could be attributed to a thermally activated expression of 
the correlation times, associated with a broad distribution of the corresponding activation 
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energies, an assumption physically plausible for many systems 

In the last twenty years, many other important contributions have advanced the under- 
standing of the 1/f noise sho wing that the presence of an 1/f spectrum can also arise from 



other basic reasons 
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4=51 ] . Thus, the con- 
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41] . is that a unique, universal origin 



elusion, now largely accepted in the literature 

m 

of 1/f does not exist, even though classes of systems can share a common basic origin of 
1/ / noise. For example, spin-glass models [l|, Q, S] provide a good explanation of the 1/f 
noise in conducting random magnetic materials. Dissipative self-organised criticality (SOC) 



models 



36 



37] clarify the origin of 1/f spectra in certain dissipative dynamical systems 



naturally evolving into a critical state. Avalanche models 



percolative models 



fl Q, Q, Q Q, y, Q 



clustering models [39] and 



represent other relevant classes of theoretical 



approaches explaining the appearence of 1/f noise in a variety of systems. In particular, 
the use of random resistor network (RRN) models III, [12|, |46j, |47|, |48] has proved to be very 



fruitful. Within this approach, a large attention has been devoted to t 



re calculation of the 



49, 



50, 



51 



52. However 



noise exponents in two-components RRN, as for example in Refs. 
these authors assumed the existence a priori of independent 1/f microscopic ffuctuators, 
associated with small fluctuations of the local resistivity, and they studied the influence 
of the topology and disorder of the network on the resistance noise magnitude at a fixed 
frequency. The existence of microscopic Lorentzian fluctuators giving rise to 1/f noise at a 
macroscopic level, has been proposed instead by Gingl et al. 53]. The hypothesis of small 
local fluctuations has been released by Seidler et al. [54] . who have shown that when the 
local resistivity fluctuations are large, the dynamical redistribution of the current gives rise 
to long-range (space) correlations and non-Gaussian 1/f noise. 

Finally, another important class of 1/f RRN models is represented by dynamical percola- 



tion models 
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41], introduced by Lust and Kakalios for describing Lorentzian spectra 



[17] and then modified to account for 1/f noise [161 . |41 |. Within these models, 1/f noise 
arises from random jumps performed by some elemental component of the system between 
two states (like trapping and detrapping of charge carriers [ly, H] , or ON and OFF states 



in the switcher model 



41] ) . Precisely, Lust and Kakalios considered a two-dimensional RRN 



with half of the resistors removed at random (percolation threshold). By focusing on the 
filamentary resistive structures connecting the electrodes, they allowed fluctuations only for 
the resistors at the nearest neighbor positions [16J, H7J]. In this manner, a uniform distribu- 
tion of the trapping times provides a Lorentzian spectrum [l^] , while a certain non uniform 
distribution of these times gives rise to an 1/f noise spectrum jl6 |. Furthermore, Celasco 



and Eggenhoffner 



41] have studied a binary RRN with links behaving as random switchers 



of resistances r, fluctuating bewteen two ON and OFF states (associated with r = and 
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r 7^ 0, respectively). The time evolution of a single switcher is controlled by two parameters, 
p and q, respectively representing the probability for a switcher to be ON at the time t and 
the probability for the same switcher to be off at t + At if it was ON at t. This model applies 
only to linear systems and it does not contain any direct link with the external conditions, 
however its nice and nearly unique feature is that it provides a resistance noise with power 
spectrum l/f a with < a < 2 depending on the values of p and q. Actually, many experi- 
ments report about a ranging between 4- 2 in the same system depending on the external 



conditions 
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231 ] . In particular, at high external biases, transitions from 1/f to 



Lorentzian noise are observed in many systems, as a consequence of the suppression of the 
plurality of characteristic times induced by extreme, far from equilibrium conditions Q, [^j] . 

The aim of this paper is to present a new model able to describe within a unified frame- 
work both Lorentzian and 1/f noises, together with the transitions between these two kinds 
of spectra 

HQ 

. Thus, by connecting the differences in the power spectral density not only 
with the parameters intrinsic to a given system but also with the different external con- 
ditions: ranging from nearly equilibrium conditions (very low biases), to non-equilibrium 
stationary states, up to failure conditions. Correspondingly, in the time domain this means 
to tune the decay of correlations from a long-term decay (power-law scaling of the two-points 
auto-correlation function) up to an exponential decay. 
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The model that we propose belongs to the dynamical percolation class [if], 0, 
Precisely, we consider a network made by several species of resistors, where each species 
is characterized by the values of the elementary resistance, of the temperature coefficient 
and of two activation energies, which control the probabilities of breaking and recovery 
processes of that species. For this reason we call this model "multi-species network" (MSN) 
model. The states, either stationary or non-stationary, of the MSN result from the stochastic 
competition between the breaking and recovery processes of the different species. In analogy 
with Dutta et al. 0|, we take the activation energies distributed in a broad range of 

values, as discussed in the next section. As a consequence, the resistance fluctuations of the 
MSN exhibit an l/f a power spectrum, where the value of a in the range -j- 2 depends on 
temperature and current. 

The paper is organised as follows: in Sect. II we illustrate the MSN model, in Sect. Ill 
we report the results and finally in Sect. IV we draw the conclusions of this study. 
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II. MODEL 



The multi-species network can 
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57|, 



58 



)e considered as a generalization of the single species 



591 ] . The SSN model describes a RN whose resistance 



network (SSN) 

fluctuates over a single time-scale. In such a network the correlations relax exponentially 
(with small deviations from a single exponential decay at high biases) and the power spectral 



density o: 
currents 



re resistance fluctuations is Lorentzian-like (a ~ 2) at all temperatures and 
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as will be discussed in the following. Both the models, MSN and 



SSN 

I 



, share some technical features that the reader can found with more details in Refs. 



57 



58]. 



As usual for RRN models 
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12 



46 



471 ]. a conducting thin film with granular structure 



is described as a two-dimensional N x N resistor network with square-lattice, where N is 
the linear size of the network [60]. The RN is biased by an external constant current, /, 
applied through perfectly conducting bars placed at the left and right hand sides and it is 
in contact with a thermal bath at a temperature T. In the MSN model proposed here the 
network is made by N spec different species of resistors, whose resistances are denoted by r n i , 
where n is the index of the elementary resistor specifying its position within the network and 
i = l,...N spec labels the species. The elementary resistances are taken linearly dependent 
on temperature: r nji = r 0i j[l + o>T,i(T n — T)], where r 0i j, ar,i are respectively the resistance 
at the thermal bath temperature and the temperature coefficient of the i-th species, while 
T n is the local temperature. By neglecting time dependent effects in the heat diffusion 61] 



we account for the local 
the local temperature as 



oule heating of the n-th resistor and its neighbors by expressing 



561 ] : T n — T + A A n , where A describes the heat coupling of the 



29]. 



elementary resistor with the thermal bath and A n depends on the local currents 

Each resistor can be in two states: regular (with resistance r n> j) or broken (with re- 
sistance tqp = 10 9 r n j). Resistors in the broken state are called defects. The transition 
from the two states is stochastic: the transition from the regular state to the broken one 
(breaking process) occurs with probability Wd,% while the reverse transition (recovery pro- 
cess) occurs with probability Wrj, where i is the species index defined above. Both pro- 
cesses are thermally activated, thus their probabilities are: Wo,i = exp(— En,i/kBT n ) and 
Wji ti = exp(—E Ri /k B T n ), where E D>i and E R)i are the activation energies of the z-species 
and ks is the Boltzmann constant. 
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The initial state of the network (equal or different fraction of each species, ordered or 
randomly distributed within the network, values of ro,j and a.T,i) has an important role on 
the network evolution: thus both on the actual achievement of stationary states and on the 
features of the resistance fluctuations. The initial conditions that we have considered in this 
work are the following. 

i) We have taken the values of ro,i and a^i uniformly distributed inside a given range of 
values: respectively r 0j j G [r mi „, r max ] and «t,j G [a m in, o: max \. This choice has been adopted 
for the sake of simplicity, while other options are reasonable as well. 

ii) We have assumed a random distribution of the species within the network. Though 
special patterns for the space distribution of the species within the network can be of interest 
in many situations, the choice adopted here is the simplest one. 

iii) We have taken the activation energies Ep,i and Er^ of the different species uniformly 
distributed inside the ranges of values [E^ m , E r ^ ax \ and [E™ n ,E^ ax ]. This choice, taken 
in analogy with Dutta et al. nag , is physically justified in many disordered systems, 
like granular and amorphous materials, composites, etc., where the orientational disorder 
present inside these materials can give rise to different energy barriers for the electron flow 



along the different conducting paths 
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12 



34|. 



vi) The energies, Eo,i, Erj have been coupled by imposing the condition that the dif- 
ference Ed^ — Erj is approximately the same for the different species: Ep^ — Er^ « AE*. 
The reasons of this assumption will be explained below. However, we anticipate that by 
defining pi as the fraction of defects belonging to the i-th species: pi = Nf, ro k,i/Ntot,i, and 7$, 
as the correlation time which characterizes its fluctuations, the conditions iii) and vi) imply 
a logarithmic distribution of the correlation times of the different species. 

v) The association between the resistance ro,j of the i-th species and the corresponding 
activation energies Ed,% and Er^ has been done by adopting the criterion that increasing 
values of r< are paired with increasing values of r ,j. Alternative choices are of course possible. 

The physical meaning and the implications of the assumptions iii) and iv) can be under- 
stood by the following arguments, which also provide a theoretical ground to the numerical 
results reported in Section III. 

Let us consider a network made by a single species of resistors at the equilibrium, in 
the vanishing current limit I — > (Ref. 1571 ). In this case, it is easy to derive the following 
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expression for the average fraction of defects under stationary fluctuations 571 ]: 



W D (1 - W R ) 1 

where Wd and are the probabilities of the breaking and recovery processes (for a SSN 
these quantities are independent of the index i and, at the equilibrium, they are also inde- 
pendent of the index n specifying the position within the network). Now, it is convenient 
to introduce the following definition: A = \o.{Wr/Wd) = (E D — E^/ksT. Since Wr <C 1, 
Eq. (0Q) can be approximated as: 

Since A controls the average fraction of defects at the equilibrium, determining the "intrinsic" 
disorder of the network, it can be called intrinsic disorder parameter. Thus, the energy 
difference Ed — Er is the effective activation energy which sets the average defect fraction 
in the network. Furthermore, we note that a SSN can be considered as a network subjected 
to random telegraph noise (RTN) of the elemental resistors, where each resistor r fluctuates 
between two states, 1 (active) and 2 (broken), with the following transition rates: W(l 
2) = W D and W{2 -»• 1) = W R /(1 - W R ). Thus, according to RTN theory 
spectral density of the elemental fluctuator is given by: 



2|, the power 



M/j " 1 + (27T/7V) 2 { } 

where p is the difference in the resistance of the two states and r r is the correlation time of 
the elemental fluctuator, given by: 

1 Wn Wn 

- = W(l - 2) + W(2 -> 1) = W D + W « = (4) 

T r (1 - Wr) (p) 

In Ref. [57| it has been shown that this expression holds also for the correlation time r of 
the fluctuations of the global network resistance R, at least when the average fraction of 
defects (p) is sufficiently far from the percolation threshold 
the last condition, it is: r = r r . On the other hand, it is well known that when a global 
quantity, like the network resistance, results from the superposition of several exponential 
relaxation processes with different correlation times, its spectral density can be written as 



46 



62]. In other terms, under 



[i 



*-cn=jf*'»-M TT ^ (5) 



where gR{r')dr ! specifies the contribution to the fluctuations of the global quantity by the 
elemental processes whose correlation times lie in the interval between r' and r' + dr'. In 
conclusion, since for a SSN at the equilibrium g^r') = 5{r' — r), the power spectral density 



Lq 



29 



581 ] has shown that 



of the network resistance fluctuations takes the Lorentzian form of 
The study of non equilibrium stationary fluctuations of a SSN 
many of the effects of an external current can be approximately described within a mean 
field-like framework, by considering average transition probabilities, (Wd) and (Wr) [58| . 
which depend on the bias through the average temperature: (T)j = T + 9R I 2 , where 
9 is the structure thermal resistance and Rq is the network resistance in the vanishing 
current limit. Of course, at high external bias, the distribution of currents and temperatures 
within the network becomes strongly non homogeneous, resulting in breaking and recovery 
probabilities, Wd = e.xp{—Er,/kBT n ) and Wr = exp(—En/kBT n ), strongly dependent on 
the position of the elemental resistor within the network. This effect, which gives rise to 



a filamented growth of the defect pattern, characteristic of biased percolation [29|, [581], also 
implies different correlations times for the fluctuations of the elemental resistors. However, 
this difference in the correlation times, is never so large to modify significantly the shape of 
the power spectral density of the network resistance fluctuations, which remains Lorentzian- 



like at all currents compatible with stationary states of the network [29_|, [58] . In other terms, 
the density function <?_r(t) in Eq. (jSJ) is different from zero only in a relatively small interval 
of t values. 

Coming back to a MSN, it is convenient to define for all the species % the quantities 
Aj = (Ez> ti — Eft^/hsT, which control the average fraction of defects (pi) at the equilibrium. 
Thus, the value of Aj determines the contribution of the z-th species to the "intrinsic" 
disorder of the network. Then, by taking: Aj ~ cost = AE* /ksT, we are assuming 
an approximately equal concentration of the different species at the equilibrium. In such 
situation, we expect that a mean field-like approach works and that the energy AE* plays 
the same role of effective activation energy as that played by the energy difference Ed — En 
in the SSN at equilibrium. In other terms, AE* controls the average defect fraction in 
the network. Furthermore, it must be noted that the conditions iii) and vi), coupled with 
Eqs. ([2]) and (jlj), imply a logarithmic distribution of the correlation times: g{T{) = 1/Vj, with 
Ti G [r m i n , T max ], where T min and r max define the time interval in which the correlation times 
are distributed. Since Wd,% and Wrj depend also on T n (i.e. on the external temperature and 
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on the local Joule heating) a logarithmic distribution of r< is obtained only at the equilibrium, 
in the vanishing current limit, when the local Joule heating is negligible, A n « 0, and all 
the resistors are at the same temperature: T n = T. Moreover, apart from this effect related 
to a non homogeneous current distribution, in the nonlinear regime there is an increase of 
the average temperature which modifies the average transition probabilities, changing T min 
and T max , according to Eq. P|. In conclusion, the values of T min and r max depend on the 
particular material (E™ n , E^ ax , E R VLn , E^ ax ) and on the external conditions (T and 7). 

The time evolution of the network is then obtained by Monte Carlo simulations which 
update the network resistance after a sweep of breaking and recovery processes, according to 
an iterative procedure detailed in Ref. |29|. The sequence of successive network configurations 
provides a resistance signal, R(t), after an appropriate calibration of the time scale. Then, 
depending on the external conditions and on the network parameters, the network either 
reaches a steady state or undergoes an irreversible electrical failure [29|, |30|, |58|, [59J]. This 



latter possibility is associated with the condition that the global average defect fraction 
(p) — YUiPi) reaches the percolation threshold, p c . Therefore, for a given network at a given 
temperature, a threshold current value, Ib, exists above which electrical breakdown occurs 



291 ] . For current I < Ib-, the steady state of the network is characterized by fluctuations of 



the defect fraction, Sp, and of the resistance, 5R, around their respective average values (p) 
and (R). 

All the results reported here concern networks of sizes 75 x 75 made by N spec = 15 resistor 
species. For the other parameters the following values have been used as representative of 
realistic cases: r min = 0.5 Q and r max = 1.5 £7, a m i n = 10~ 4 K™ 1 , a max = 1CT 1 K _1 . 
Moreover, we have taken: E^ in = 58 meV, E% ax = 375 meV, E% in = 37 meV, E% ax = 346 
meV and AE* ~ 25 meV (precisely, the average value of the difference between Eo^/ks 
and E Rti /k B is AE*/k B w 319.86 K). 

According to Eq. (j4j) and conditions iii) and iv), in equilibrium (or nearly equilibrium) 
at the reference temperature T re f = 300 K, the above energy values imply (p^ w 0.25 Vi, 
Tmin ~ 2 and T max si 5 x 10 5 (where times are expressed in units of iterative steps). At 
T > T re f, the interval defined by T min and r max becomes progressively narrower while (p^ 
increases. Of course, the contrary happens for T < T re f. We underline that these values 
of (p^, r min and r max , calculated from Eq. (j3J) are reported here just to give a qualitative 
idea of the network state and because they can be useful for the choice of the parameters 
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gmin ^ gmax ^ gmtn ^ grnax However all the network properties, including the average defect 
fraction, the correlation time of the resistance fluctuations and the other results discussed in 
Sec. Ill are obtained directly from the output of simulations. Finally, the auto-correlation 
functions and the power spectral densities of the resistance fluctuations are calculated by 
analyzing stationary R(t) signals consisting of 1 -f- 2 x 10 6 records. 

III. RESULTS 

Figure 1 reports the resistance evolution of a MSN calculated at 300 K in the vanishing 
current limit. The inset displays a small part of the same evolution over an enlarged time 
scale. Here we notice the co-existence of different characteristic time scales in the R(t) signal. 
Indeed, the long relaxation time associated with the achievement of the steady state, r re /, 
co-exists with the shorter times characterizing the resistance fluctuations, as displayed in the 
inset. For comparison, Fig. 2 reports the time evolution of the resistance of a SSN obtained 
at T = 300 K in the same bias conditions. In this case, the values of the activation energies 
(E D = 350 meV and E R = 310 meV) are chosen to give a relaxation time comparable with 
that of the signal in Fig. 1. Now, the resistance signal is controlled by a single time scale 
(r « T re i) and it is essentially flat over time scales shorter than r. This is emphasized by 
the inset in Fig. 2 where the stochastic signal resembles that of a few levels system (it must 
be noted that the vertical scale of the inset in Fig. 2 is significantly enhanced with respect 
to that of Fig. 1). 

Figure 3 displays the auto-correlation functions of the resistance fluctuations in Figs. 1 
and 2, corresponding to the MSN model (black triangles) and to the SSN model (black 
squares). A log-log representation is adopted for convenience. The solid and short dashed 
lines represent the best-fits to the two auto-correlation functions carried out, respectively, 
with a power-law and an exponential law. The fitting procedure confirms the exponential 
decay of the correlations in the resistance fluctuations of the SSN. Furthermore, it points 
out the existence of long-term correlations in the resistance fluctuations of the MSN, char- 
acterized by a power-law decay of the auto-correlation function: 

C SR (t) ~ (6) 

with < 7 < 1. In particular, here we have found a value 7 = 0.22 ±0.01 for the correlation 
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exponent. We stress that the above expression for the auto-correlation function implies a 
divergence of the correlation time, as can be easily seen by considering the following general 
definition of the correlation time 64|: 



CsR ® dt (7) 

Figure 4 shows the power spectral densities of the resistance fluctuations of a MSN and of 
a SSN calculated at 300 K by Fourier trasforming the auto-correlations functions of Fig. 3. 
Here, the two grey solid lines represent the best-fits with a power-law of the MSN spectrum 
in the low and high frequency regions. We can see that the resistance fluctuations of the 
multi-species network exhibit at low frequencies a power spectral density scaling as l/f a , 
with a value a = 0.94. We notice that this scaling behavior holds over several decades of 
frequency. This result is a consequence of the envelope of the different time scales associated 
with the different resistor species, described by Eq. (j5J). Instead, in the high frequency region, 
the slope of the spectrum is -1.53 (in fact, the slow relaxations are ineffective at such high 
frequencies). For contrast, the grey dashed curve in Fig. 4 is the best-fit with a Lorentzian 
to the SSN spectrum. The corner frequency of the Lorentzian, f c = 4.0 x 10~ 6 (arbitrary 
units), is consistent with the correlation time reported in Fig. 3 and obtained by the best-fit 
of the corresponding auto-correlation function. 

Now, we will discuss how the temperature of the thermal bath affects the resistance 
fluctuations of a multi-species network at equilibrium or nearly equilibrium conditions (low 
biases). In other terms, we will analyse the properties of the fluctuations in the Ohmic 
regime and for different temperatures. Figure 5 reports the resistance evolutions calculated 
at increasing temperatures: T = 400 K (lower curve) and T = 600 K (upper curve). Already 
this qualitative comparison between the R(t) signals shows that a temperature increase 
implies a significant growth of both the average resistance and the variance of the resistance 
fluctuations (see also Figs. 1). Furthermore, this comparison points out a drastic reduction 
of the relaxation time at increasing temperatures (more than one order of magnitude when 
T rises from 300 K to 400 K). 

The temperature is also found to affect the distribution of the resistance fluctuations by 
increasing its skewness, as shown in Fig. 6, which reports the probability density function 
(PDF) of the resistance fluctuations for several temperatures. A normalized lin-log represen- 
tation has been adopted here for convenience (a is the root mean square deviation from the 

11 



average resistance). The normalized Gaussian distribution is also reported for comparison. 
The figure shows a significant non-Gaussianity of the resistance fluctuations, which becomes 
stronger at high temperatures, when the system approaches failure conditions. This behav- 
ior is completely different from the behavior of the defect fraction fluctuations, shown in 
Fig. 7, which remains Gaussian at all temperatures. This different behavior can be easily 
understood. Actually, at increasing temperature, the local resistivity fluctuations become 
progressively larger. This fact emphasizes the dynamical redistribution of currents and 
brings to the emergence of long-range correlations inside the network 54( . Then, the viola- 



tion o 
noise 



' the vali dity conditions of the central limit theorem leads to a non-Gaussian resistance 



19. 



54 



651 ] . In other terms, the results in Fig. 6 imply that the correlation length 



of the resistance fluctuations progressively increases with the temperature. A discussion of 
the distribution of the resistance fluctuations and the role played on the non-Gaussianity 
by the size, shape and disorder of a SSN is reported in Refs. 



30, 



59 



where the link with 



;he universal distribution of the fluctuations of Bramwell, Holdsworth and Pinton (BHP) 



66 



671 ] is analyzed. 



Figure 8 reports the auto-correlation functions of the resistance fluctuations of a MSN 
calculated at 400 K and 600 K. To help the comparison, the auto-correlation function at 300 
K (already shown in Fig. 3) has been drawn again in Fig. 8, together with its power-law 
best-fit (grey solid line). The dashed grey curves represent the best-fit to the correlation 
functions at 400 K and 600 K with the expression: 



C SR (t)=C t- h exp(-t/u) 



(8) 



The values of the best-fit parameters are: Co = 1.13, h = 0.30 and u = 1.42 x 10 4 for 
T = 400 K and C = 0.965, b! = 0.46 and u' = 5.57 x 10 2 for T = 600 K. The fit to 
Csr with Eq. ([8]) is found to be very satisfactory. We conclude that at T > 300 K, the 
auto-correlation function of the resistance fluctuations of the MSN is well described by a 
power-law with an exponential cut-off. Such kinds of mixed decays of the correlations, non- 
exponential and non power-law, are often found in the transition of a complex system from 
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55 



68 



69, 



70 



3, 0,0- it 



short-term correlated to long-term correlated regimes 
should be noted that for u — > oo, Eq. ([8]) becomes a power-law, while for h — > 0, it describes 
an exponential decay. Actually, Fig. 8 highlights a strong reduction of the correlation time 
of the resistance fluctuations as the temperature increases. This result can be understood 
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in terms of Eq. (j4]) , by considering the thermally activated expressions of the breaking and 
recovery probabilities. Indeed, the increase of temperature above T re f implies the reduction 
of the ratio r min /T max and the correspondent narrowing of the interval [T min ,T max ], where 
the Tj are distributed. This trend tends to suppress the power-law decay of correlations in 
favor of the exponential decay. For similar reasons, the temperature decrease below T re f, 
implies the opposite trend with the correlations keeping their power-law decay over wider 
time scales. For example, at T = 200 K the ratio r m j n /r max becomes 7.5 x 10 7 . Thus, we 
can define the transition temperature, T*, as the temperature value which signs the crossing 
from a long-term correlated behavior (occurring for T < T*), to a behavior characterized by 
a finite and relatively short correlation time (occurring for T > T*). 

The correlation time of the resistance fluctuations can be directly estimated by making 
use of Eq. ([7j) and Eq. ([8]). In terms of the best-fit parameters of the auto-correlation 
function, it is easy to derive the following analytical expression for r: 

t = u^Til - h) (9) 

where V is the Gamma function. The values of r calculated in this way increase at decreasing 
temperatures. In particular, r exhibits a sharp increase when T approaches T*, in agreement 
with the long-term decay of the correlation function at this temperature. We have found 
that the behavior of r versus temperature is well fitted by the power-law: r ~ (T — T*)~ e . 
Then, the fit procedure allows us to determine the value of the transition temperature. We 
have found: T* = 306 K and 9 = 2.7. Figure 9 reports the values of r as a function of the 
difference T — T*. The dashed straight line corresponds to the above mentioned power- law. 
Therefore, we can conclude that: T* w AE* /ks- Furthermore, as a result of the particular 
choice of the parameters adopted in the present calculations, the transition temperature is 
close also to the reference temperature, T* « T re f. However we remark that the temperature 
T re f has been introduced in the model merely to help the choice of the activation energies (to 
provide a sufficiently wide interval [r min , T max ]). While, AE* is the relevant input parameter 
which determines T*. 

Figure 10 displays the spectral densities of the resistance fluctuations calculated at 400 
K and 600 K. The grey lines represent the best-fits with a power-law of the MSN spectra 
in the low and high frequency regions. Precisely, at low-frequencies, the slopes of the lines 
are -0.87 and -0.78 respectively for T = 400 K and T = 600 K. Instead, at high frequencies 
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the respective slopes are -1.31 and -1.02. Thus, in both frequency regions, low and high, 
the slopes of the specra are reduced when the temperature increases. We conclude that at 
T > T* the power spectrum keeps the 1/ f a form with the value of a significantly decreasing 
below unity. This decrease is pointed out in Fig. 11, which reports a as a function of the 
temperature. The dashed line in Fig. 11 is the best-fit with a linear law. We notice that 
such a decrease of the exponent, from a ~ 1 to a ~ 0.8 — 0.5, is frequently observed in 
the experiments at intermediate temperatures On the other hand, many 

experiments have pointed out a strong dependence of the detailed behavior of a(T) also 



on the particular material 



flfl 



181 ]. In this respect, we remark that in this work we 



are focusing our attention to the general features of the spectra, related to the term of 
correlations, rather than to the interpretation of a particular set of experiments. Actually, 
in our model the quantity AE* (determining the transition temperature T*) is the input 
parameter which can be adjusted for a quantitative fit of experiments. Furthermore, we 
stress that the monotonic decrease of a for T > T* is obtained here in the linear regime of 
currents, i.e. neglecting Joule heating effects. Actually, these effects, whose importance also 
depends on the temperature, can give rise to a more complicated behavior of a versus T, 
which are outside the interest of this research but can be implemented with minor efforts. 

Now, to complete the analysis of the MSN model, we consider the dependence on the 
temperature of the average resistance, (R), and of the variance of the resistance fluctuations, 
((AR) 2 } for a multi -species network in nearly equilibrium conditions. This dependence has 
been already qualitatively depicted in Figs. 1 and 5. However here we want to quantify 
these behaviors, also with the purpose of checking by numerical results the validity of the 
mean field-like approach discussed in Section II. Figure 12 displays (R) as a function of the 
temperature, while the inset shows the dependence on the temperature of the average defect 
fraction. We have performed a one parameter best-fit of the numerical values of (p) with the 
expression: (p) = c p /(l + 6 ae */ k bT^ ^y e h ave f oun d c p = 0.93 ± 0.03, a value compatible 
with the expression: (p) « 1/(1 + e AE */ K BTy Thus, Eq. (TSJ) qualitatively accounts for the 
behavior of the average defect fraction versus temperature. 

For what concerns the dependence on temperature of the average resistance, in percolation 
theory 46( it is well known the following scaling relation between the network resistance 
and the defect fraction: R ~ \p — p c | _At . We note that in the present model at the vanishing 
current limit j^lj], the percolation is uncorrelated and the exponent \i takes the universal 
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value [46| \i = 1.303, while p c = 0.5 [46j, |62(. The dashed curve in Fig. 12 shows the best-fit 
of the numerical data for the average resistance with the expression: (R) = cr\ (p) — p c \ In 
the best-fit procedure the values of and p c have been taken fixed to their theoretical values, 
while we have taken: (p) = c p /(l + e AE */ K BTy va ] ues found for the fitting parameters 
are: cr = 0.42 ± 0.01 Q and c p = 0.80 ± 0.03. Thus, the dependence on temperature of 
the average resistance of a MSN in the linear regime is well described by the usual scaling 
relation coupled with the mean field-like expression of (p). 

Figure 13 reports the variance of the resistance fluctuations as a function of the temper- 
ature. This behavior can be easily understood by considering before the power-law rela- 



tion between the variance of the resistance fluctuations and the average resistance 46j, |49) : 
((AR) 2 ) ~ (R) v , often written in terms of the relative variance, as: 

((AR) 2 )/(R) 2 ~ (R) s (10) 

where s = rj — 2. By reporting on a log-log plot the relative variance ((AR) 2 ) / (R) 2 versus 
(R), the calculated slope provides for the exponent s the value s = 2.6 ± 0.01 We note 



that this value of the relative noise exponent agrees with the value reported in Ref. 



57 



This agreement is consistent with the fact that the MSN model generalizes the results of 



Ref. l57j to networks characterized by l/f a noise. Now, we can use the information on the 
exponent 77 = s + 2 = 4.6 to check the consistency among the dependence on temperature 
of ((AR) 2 ), Eq. ([2]), the scaling relations between (R) ~ (p) and between ((AR) 2 ) ~ (R). 
The solid curve in Fig. 13 shows the best-fit to the variance of the resistance fluctuations 
considered versus temperature, with the following expression: ((AR) 2 ) = c^c^Kp) — p c 



where the fitting parameters take the values: ca = (4.42 ± 0.04) x 10 , cr = 1.0 ± 0.04 
and c p = 0.74 ± 0.03. Thus, a mean field-like framework is overall able to account for the 
dependence on temperature of the lowest two moments of resistance fluctuation distribution. 

Until now we have analysed the resistance fluctuations of the network in equilibrium 
or in nearly equilibrium conditions, in the Ohmic regime, and for different temperatures. 
Now, we will briefly discuss the non equilibrium properties of the fluctuations at high biases, 
in the nonlinear regime, and for a given temperature. Figure 15 shows the auto-correlation 
functions of the resistance fluctuations of a MSN for increasing values of the external current. 
All the functions are calculated at 300 K. The open squares represent the auto-correlation 
of the resistance fluctuations in the Ohmic regime: 1 = 5 mA, the other three curves 
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correspond to the non-linear regime. Precisely, black triangles: I = 200 mA; open circles: 
/ = 220 mA; grey diamonds: I = 250 mA > Is- Thus, the last curve corresponds to failure 
of the network, i.e. non-stationary resistance fluctuations. In this case, the auto-correlation 
function has been calculated by considering only the nearly-stationary portion of the R(t) 
signal (after subtraction of a linear trend). The dotted grey line is the best-fit with a power- 
law of slope —0.22 ± 0.01, the dashed lines are the best-fits with Eq. ([8]). The values of 
the parameter h corresponding respectively to: I = 200, 220, 250 mA are the following: 
h = 0.19, 0.17, 0.11 (in all cases, the numerical error is estimated as ±0.01). The values of 
the parameter u corresponding to the same currents respectively are: (1.26 ± 0.01) x 10 5 , 
(9.15 ± 0.01) x 10 4 , (5.84 ± 0.01) x 10 3 . Thus, both h and u decrease at increasing biases 
and for high biases h — > 0. In other terms, at increasing currents, time correlations decay 
progressively faster, the long-term correlated behavior is destroyed and it emerges a trend 
towards a simple exponential decay. 

The corresponding effect in the frequency domain of increasing biases is shown in Fig. 16, 
which reports the power spectral densities at 300 K. The curve (1) is obtained for the Ohmic 
regime (1 = 5 mA); the other three curves display the spectra in the non-linear regime. The 
solid grey lines are the best-fits with power-laws in the low and high frequency regions of 
the spectrum. The slopes are specified in the figure. We note that: i) the corner frequency 
between the two regions progressively moves towards lower frequencies at increasing biases; 
ii) the slopes in the high frequency region are: —1.56, —1.62, —1.66, —1.78 respectively for 
1 = 5, 200, 220 and 250 mA, i.e. the slopes increase at increasing currents (in all cases, the 
numerical error on the slope values is estimated as ±0.01). Thus, at high biases, the spectral 
densities show a trend towards a transition from an 1// to a Lorentzian behavior. However, 
we remark that for the present choice of the initial conditions and/or the numerical values 
of the parameters used in these calculations we have not got a pure Lorentzian spectrum. 

IV. CONCLUSIONS 

We have developed a stochastic model to investigate the l/f a , with < a < 2, re- 
sistance noise in disordered materials. More precisely, we have considered the resistance 
fluctuations of a thin resistor with granular structure in different stationary states: from 
nearly equilibrium up to far from equilibrium conditions. Furthermore we have also consid- 
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ered fluctuations in non-stationary states, associated with failure of the electric properties. 
This system has been modeled as a two-dimensional network made by different species of 
elementary resistors. The steady state of this multi-species network is determined by the 
competition among different thermally activated and stochastic processes of breaking and 
recovery of the elementary resistors. The network properties have been studied by Monte 
Carlo simulations as a function of the temperature and applied current, in both Ohmic and 
non-Ohmic regimes. A mean field-like framework has been also used to qualitatively describe 
the dependence on temperature of the lowest two moments of the resistance fluctuation dis- 
tribution. Furthermore, the correlation properties of the resistance fluctuations have been 
analyzed in both the time and the frequency domains. The model gives rise to resistance 
fluctuations with different power spectra, depending on the external conditions. Thus it 
provides a unified approach to the study of materials exhibiting either Lorentzian noise or 
1/ f a noise. By analyzing the correlations in the time domain, it has been found that the 
resistance fluctuations display a cross-over from long-term correlations to intermediate-term 
correlations. Although a trend towards an exponential decay has been identified, for the 
present choice of the parameters/initial conditions, we have not got a single exponential de- 
cay. However, the model proposed seems able to account for the complex interplay exerted 
on the correlation properties of the resistance fluctuations by the external conditions and 
resulting in the complicate behavior of the noise exponent observed in many experiments. 
In perspective, it should be explored the role of different initial conditions: such as unequal 
presence of the species and/or non-uniform distribution inside the network, which can give 
rise to interesting spatio-temporal organization patterns. 
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FIG. 1: Resistance evolution of a multi-species network (MSN model) calculated at 300 K. The 
resistance is espressed in Ohm and the time in iterative steps. The inset highlights the resistance 
fluctuations on an enlarged time scale. In particular, the time units are divided by a factor 10~ 5 . 
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FIG. 2: Resistance evolution of a single-species network (SSN model) calculated at 300 K. The 
resistance is espressed in Ohm and the time in iterative steps. The inset displays the resistance 
fluctuations on the same enlarged time scale of the inset in Fig. 1. (the vertical scales of the insets 
in the two figures are different). 
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FIG. 3: Auto-correlation functions of the resistance fluctuations calculated for a MSN (black 
triangles) and for a SSN (black squares). Both functions are obtained at 300 K. The solid and 
short dashed grey lines show the best-fit respectively with a power-law of exponent 7 = 0.22 and 
with an exponential with correlation time r = 5.8 x 10 4 . The time is expressed in iterative steps. 




FIG. 4: Power spectral density of the resistance fluctuations at 300 K calculated for a MSN (solid 
line) and for a SSN (dotted line). The grey solid line shows the best- fit to the MSN spectrum with 
a power-law of slope -0.94. The grey dashed curve represents the best-fit with a Lorentzian to the 
SSN spectrum. 
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FIG. 5: Resistance evolution of a MSN at 400 K and 600 K. The resistance is espressed in Ohm 
and the time in iterative steps. 




(<R>-R)/o 

FIG. 6: Normalized probability densities of the resistance fluctuations of a MSN calculated at 300 
K (full diamonds), 400 K (open circles), 500 K (open up-triangles) , 600 K (full down-triangles) and 
700 K (stars), a is the root-mean-square deviation from the average resistance. The solid black 
curve is the Gaussian distribution and the dashed one the BHP distribution (see text). 
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FIG. 7: Normalized probability densities of the defect fraction fluctuations calculated at 300 K (full 
diamonds), 400 K (open circles), 500 K (open up-triangles) , 600 K (full down-triangles) and 700 
K (stars). Here a is the root-mean-square deviation from the average value of the defect fraction. 
The solid black curve is the Gaussian distribution. 
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FIG. 8: Auto-correlation functions of the resistance fluctuations of a MSN calculated at different 
temperatures. The solid grey curve shows the best-fit with a power-law to the auto-correlation 
function at 300 K (the same of Fig. 3). The dashed grey curves display the best-fit to the auto- 
correlation functions at 400 and 600 K with the function: C(t) = Cot~ h exp[— t/u] (see the text for 
the values of the fit parameters). The time is expressed in iterative steps. 
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FIG. 9: Correlation time of the resistance fluctuations of a MSN as a function of the difference 
T — T*. The time is expressed in iterative steps and the temperature in K. The value of T* is 
reported in the figure. The dashed line shows the fit with a power-law of exponent 9 = 2.7. 
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FIG. 10: Power spectral density of the resistance fluctuations of a MSN at T = 400 K and T = 600 
K. The grey lines show the best-fit with power-laws of slopes -0.87 and -0.78, respectively. The 
spectral density at 600 K has been multiplied by a factor 10 for visual reasons. 
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FIG. 11: Noise exponent a as function of the temperature. The dashed line refers to a linear 
best-fit of the calculated values of the exponent. 
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FIG. 12: Variation of the average resistance (full circles) and average defect fraction (stars, in the 
inset) as function of the temperature. The resistance is espressed in Ohm and the temperature in 
K. The solid and dashed straight lines are the best-fits with the expressions reported in the figure 
(see the text for the values of fit parameters). 
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FIG. 13: Variance of the resistance fluctuations as function of the temperature. The variance 
is espressed in and the temperature in K. The dashed straight line shows a best-fit with the 
expression reported in the figure (see the text for the values of fit parameters). 
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FIG. 14: Relative variance of the resistance fluctuations as a function of the average resistance (this 
last expressed in Ohm). The dashed line shows a best-fit with a power-law of exponent s = 2.6. 
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FIG. 15: Auto-correlation functions of the resistance fluctuations of a MSN at room temperature 
for increasing value of the external current. Ohmic regime: 1 = 5 mA (open square), non-linear- 
regime: / = 200 mA (black triangles), / = 220 mA (open circles) and / = 250 mA (grey diamonds). 
The dotted grey line is the best-fit with a power-law, the dashed lines are the best-fit with the 
function: C(t) = Cot~ h exp[— t/u] (see the text for the values of the fit parameters). 



FIG. 16: Power spectral density of the resistance fluctuations of a MSN at room temperature for 
increasing value of the external current. Ohmic regime: curve (1) / = 5 mA; non-linear-regime: 
curve (2) / = 200 mA, curve (3) I = 220 mA and curve (4) / = 250 mA. The solid grey lines 
are the best-fit with power-laws, the resulting slopes in the different regions of the spectrum are 
reported in the figure. For visual reasons curve (2) has been multiplied by a factor 2 x 10 2 , curve 
(3) by a factor 2 x 10 4 and curve (4) by 1 x 10 6 . 
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